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ABSTRACT 

We derive distance, density and metallicity distribution of the stellar Monoceros Overdensity (MO) 
in the outer Milky Way, based on deep imaging with the Subaru Telescope. We applied CMD fit- 
ting techniques in three stripes at galactic longitudes: / ~ 130°, 150°, 170°; and galactic latitudes: 
+15° <6<+25°. 

The MO appears as a wall of stars at a heliocentric distance of ~10.1±0.5 kpc across the observed 
longitude range with no distance change. The MO stars are more metal rich ([Fe/H]~ —1.0) than the 
nearby stars at the same latitude. These data are used to test three different models for the origin 
of the MO: a perturbed disc model, which predicts a significant drop in density adjacent to the MO 
that is not seen; a basic flared disc model, which can give a good match to the density profile but the 
MO metallicity implies the disc is too metal rich to source the MO stars; and a tidal stream model, 
which, from the literature, bracket the distances and densities we derive for the MO, suggesting that 
a model can be found that would fully fit the MO data. Further data and modeling will be required 
to confirm or rule out the MO feature as a stream or as a flaring of the disc. 
Subject headings: Galaxy:disk,Galaxy:structurc 



1. INTRODUCTION 

The Monoceros Overdensity (MO) is an extensive stel- 
lar structure found in the outer regions of the Milky Way 
at Galactocentric distances of ~15 - 18 kpc. It was origi- 
nal ly discovered in t he Slo an Digital Sky Survey (SDSS) 
by Newberg et al. (2002) and subsequent observations 
reveal a similar structure in many directions around the 
Galaxy ( Yanny et al. |2003[ Ibata et al. 2003| Crane et 



aT^[2Ti03linonn et al.||2005a|bl |Martln et ■ al.|2Wtia||Conn 
et al.||20071|2008[ |Casetti-Dinescu et al.||2008[|Solhma et 
l.| |201ip ~From this, it has been concluded that it forms 



a coherent structure from Galactic longitudes of I = 60° 
- 300° and straddles both sides of the Galactic plane. 
While the approximate extent of the MO is tentatively 
mapped, its origins are somewhat obscure. 
The MO formation scenarios fall into three broad cate- 



gories: tidal t ails from an accreting dw arf galaxy ( Martin 



et al. 2004a Peharrubia et al. 20051; misidentification 



ot normal Galactic warp /Hare profiles (Momany et al 



20041 |Moitinho et al.|[2006l |Momany et al.||2TO| |L6pez- 
Corredoira et al.||2007[ IHammersley fc Lopez-Corredoira 



201ip ; and perturbed disc scenarios whereby a close 



encounter with a massive satel lite induces rings in the 



outer disc from local material (Kazantzidis et al. 



2008 



|Michel-Dansac et al.|[20TT] |G6mez et al.||iMi| |Purcell et 
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we consider each scenario briefly then the first sce- 
nario is a Galactic accretion event, where the MO is en- 
visioned to be the tidal tails of a dwarf galaxy merging 
in-Plane with the disc. Such a scenario is attractive since 
it links with the A-Cold-Dark-Matter cosmol ogy where 



galaxies for m via successive accretion events (White & 



Rees 1978) and the discover y of many stellar streams 



in the halo of the Milky Way jNewberg et al.1|2002l IbT| 



lokurov et aL"1|2006| |Grillmair pM6b[ |Ornlmair[|26b9) . 



Lhe proposed progenitor for the MO is the putative Ca- 



ms Ma jor dwarf galaxy, first discussed in [Martin et al 
( 2004a ) . If t rue , the MO could be a relic of formation as 



discussed in Freeman fc Bland-Hawthorn | ( |2002[ ) . 

The second scenario is used to explain both the stel- 
lar overdensity in Canis Major and the MO in terms of 
standard properties of a galactic disc, that is, the warp 
and the flare. The Milky Way disc is observed, with var- 
ious tracers, to warp up in the first two quadrants and 
warp down in the second two quadrants. For instance, 



2002 ]) follows th is feature with 



Lopez-Corredoira et al. 

red clump giant stars while Yusitov ( 2004 ) uses pulsars as 
tracers. As the density of the disc drops with increasing 
galactic radii, it thickens and flares. The Canis Major 
dwarf galaxy is therefore the disc of the Galaxy dipping 
below the plane and the MO is the intersection of the 
flaring disc at high latitudes above the plane. 

Finally, the third scenario invokes the interaction be- 
tween a dark matter satellite and disc to induce the for- 
mation of rings at large galactic radii. The repeated 
passage of these satellites through the disc drives the 
formation of spiral arms and rings. This has been tested 



for lo w inclination satellite encounters ( Kazantzidis et al. 
|2008|) and more recently for the Sagittarius dwart galaxy 



ichcl -Dansac et al.|2011 Gomez et al.|2011 Purcell et 



aT 201l[ ), which is on a polar orbit. 



All ot these scenarios are able, with varying levels of 
success, to fit the general spatial and kinematic profiles 
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of the MO and previous attempts to find some decisive 
evidence or prediction has not ruled out any of these 
possibilities. More and more observations of the MO, 
both photometric and spectroscopic are becoming avail- 
able (SEGUlQ PanSTARRS^l SkyMapper, etc) and so 
these degeneracies may be broken soon. 

In order to shed light into the above dilemma, we se- 
cured observations using the Subaru telescope. Our goal 
is to investigate the spatial density profile of the MO and 
also to compare our observational results with the differ- 
ent theoretical scenarios presented in the literature. In 
Section [2] we discuss the data preparation in terms of the 
observations, reduction and calibration of the dataset. 
Section [3] outlines the analysis of the color magnitude 
diagrams using CMD fitting techniques. In Section H] 
we compare our findings with the current scenarios of 
formation for the MO and in Section [5] we present our 
conclusions. 

2. DATA PREPARATION 

2.1. Observations and Reduction 

The data were collected using the SUPRIME-CAM 
Wide Field Imager on the Subaru Telescope in Hawaii. 
SUPRIME-CAM is a 10 chip camera with a field of view 
of 34' x 27' and a pixel scale of 0".20. These observa- 
tions took place in Service mode and were carried out on 
the 9 th of November 2007 and the 1 st of January 2008. 
These observations are summarized in Table [l] Roughly, 
180 frames were taken in two filters, Sloan g and r, and 
arranged into the 3 stripes across the thick disc of the 
Galaxy. Each g frame was 124 seconds and each r frame 
was 76 seconds. Figure [T] shows the survey layout for 
each stripe with the location of the fields depicted as red 
polygons overplotted on the local dust extinction con- 
tours. Implementing the program in this manner with 
SUPRIME-CAM, allows for deep observations to be ob- 
tained across large areas in a short period of time as 
required by this study. The final survey locations are 
the result of observational constraints and data quality 
issues. 

The data were reduced using the Cambridge Astro- 
nomical Survey Unit Wide Field Camera Pipeline ( Irwin 
fe; Lewis p001[ ). This pipeline was originally developed 
for the Isaac Newton Telescope Wide Field Camera and 
has since been modified to work on most of the Wide 
Field Imagers available today. The pipeline reduced data 
has been bias-subtracted, flat fielded using twilight sky 
flats and then flat fielded again using a dark sky super- 
flat. Following this, the photometry and astrometry have 
been determined using the same pipeline. The accuracy 
and co mple teness of the photometry will be discussed in 
Section 2.4 The astrometry, based on the 2MASS point 
source catalogue, is typically accurate to between 0.2 and 
0.3 arc seconds. 

2.2. Correcting the Photometry using SDSS 

Having generated the catalogue of sources and classi- 
fying them, the photometric calibration was performed 
by cross-matching sources with t he SDSS Data Release 6 
( jAdelman-McCarthy et al.|2008| . At the time of the sur- 
vey only a few of the fields overlapped with the SDSS and 

8 Sloan Extension for Galactic Understanding and Exploration 

9 Panoramic Survey Telescope And Rapid Response System 



TABLE 1 
Summary of Observations. 



Region 


Fields 


Date Observed 


Filter 


Median Seeing 


130 stripe 


1 - 16 


2007-11-09 


9,r 


0.71" ,0.62" 




20 - 32 


2008-01-08 


9 


0.78" 




17 - 29 


2008-01-08 


r 


0.69" 


150 stripe 


1 - 32 


2007-11-09 


9 


0.59" 




1 - 32 


2008-01-08 


r 


0.74" 


170 stripe 


1 - 32 


2007-11-09 


9,r 


0.47" ,0.5" 



so the calibration was performed on those fields and then 
applied to the others according to their date of observa- 
tion. Two offsets were needed as it was noticed that Chip 
10 has a much lower efficiency than the rest of the array. 
Since the 150 stripe also had SDSS corrections available, 
when correcting the other fields the choice between these 
two was based on which night those observations were 
taken on. The final photometric solution has a typical 
scatter of 3.5% in the g-band and 2.4% in the r-band 
around the SDSS values. 

2.3. Correcting the Photometry using dust extinction 

maps 

After the initial correct ion using SDSS, a correction 
b ased o n the dust maps of Schlegel, Finkbeiner, fc Davis 
| ( 19981) c ombined with the adjustment of Bonifacio et 
ah 12000) was implemented. The final photometry is pre 
sented as Hess CMDs in Figure [2] In general, the dust 
contamination is less than an E(B — V) of 0.2, for the 
majority of the survey. The data here has been extinc- 
tion corrected and represents all the fields for each stripe 
combined into one figure. 



2.4. Magnitude Completeness 

For those fields which overlap with each other, the com- 
pleteness of this sample has been determined in the same 
manner as used in the 2MASS All-Sky Point Source Cat- 



alogue (Skrutskie et al. 



2006 1 . By determining the frac- 



tion of stars that are detected in both overlapping im- 
ages as a function of magnitude with respect to the total 
number of stars observed, an estimate of the complete- 
ness can be made. This photometric completeness curve 
is then fit by the Logistic function: 



CF = 



1 



(1 + e( m_ ™')/ A ) 



(1) 



where m is the magnitude of the star, m c is the magni- 
tude at 50% completeness and A characterizes the width 
of the rollover from 100% to 0% completeness. The aver- 
age values found in each field are presented in Table [2| a nd 
an example completeness profile is shown in Figure |3J 

2.5. Survey Field CMDs 

The pipeline provides a set of classifications for dif- 
ferent objects. We are interested in those classified as 
stellar and possibly stellar object s. The detailed informa - 
tion on the categories is given in Irwin & Lewis (2001 ). 
In short, each processed frame m the pipeline is ana- 



lyzed using an object detection algorithm based on Irwin 



( |1985||1997[ ). Generated parameters include information 
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Fig. 1. — Survey layout against Galactic Extinction from the dust maps from Schlcgcl, Finkbeincr, & Davis (1998). The red polygons 
represent the position of the fields surveyed here. The final survey locations are the result ot observational constraints and data quality 
issues. The grey scale shows the underlying dust extinction with the blue contours outlining the E(B — V) values of 0.05, 0.1, 0.2, 0.4 and 
0.6. 



TABLE 2 
Mean Magnitude at the 50% 

COMPLETENESS LEVEL (m c ) AND 
TYPICAL WIDTH OF THE ROLLOVER 
FUNCTION (A) FOR EACH FILTER. SEE 

Equation [l] 



Region 


m c (g a ) 


m c (r D ) 


A 


130 Stripe 


24.2 


23.2 


0.73 


150 Stripe 


24.8 


23.8 


0.55 


170 Stripe 


25.0 


24.2 


0.70 



on position, intensity and shape. To discriminate be- 
tween background galaxies and real objects three flux 
estimates are made: integration of the flux above the 
specific age; the detection isophote for each image is ex- 
panded using an elliptical aperture to perform a curve- 
of-growth analysis; and a 'poor man's' PSF fit using a 
radius equal to the FWHM. The stellar objects and pos- 
sibly stellar objects are then selected from all the fields 
within a given stripe and plotted in a single CMD. Fig- 
ure [2] shows the three (g _r o,<?o) deep Hess CMDs with 
each containing approximately 200,000 stars which cor- 
respond to the three longitudinal stripes of the survey. 
To allow the low density MO feature to be clearly visible, 
the CMD greyscale is scaled using the square root of the 



counts shown in each pixel. These exceptional CMDs, 
reaching more than three magnitudes below the oldest 
MS turnoffs (g ~ 19.5), are the deepest observations of 
the MO to date. The CMDs show the old MS population 
along its complete extent from the blue turn-off region to 
faint red MS stars. The high quality of the photometry 
and the small errors on the main sequence allow us to 
secure the distance-metallicity degeneracy in the CMDs. 
This ensures we can quantitatively measure the stellar 
content of the MO. 

3. DATA ANALYSIS 

3.1. CMD Fitting Technique 

In order to obtain the stellar populations' structure 
at the location of the Monoceros Overdensity w e used 
the MATCH software package polphin| [2002] ) in its 
distance-fitting mode. MATCH was originally developed 
to obtain quantitative star formation histories (SFHs) 
and age-metallicity relations for systems in which all the 
stars are assumed to be at the same distance. For this 
purpose, the distance is fixed and the age and metallicity 
are independent variables. In this paper, we apply the 
CMD-fitting techniques to span the local stellar popu- 
lations within the Milky Way. For this goal, we can no 
longer consider that all stars are equidistant and, thus, 
the distance is a free parameter. In order to limit the 
number of free parameters, we define a set of template 
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1 = 130° 1 = 150° 1 = 170° 




Fig. 2, — Deep Hess CMDs of p oint sources obtained for each of the three stripes. The data has been extinction corrected according 
to |Schlegel, Finkbeiner, fe Davis"] | |1998| l and |Bonifacio et al.| | |2000) . The cloud of sources in the bottom-left of each panel stems from 
unresolved galaxies. A difference in the image quality between the observations of I = 130° and the I = 170° stripe is responsible for the 



increase in depth of the resulting CMDs (see Table Main sequence turn-off (MSTO) stars, which we associated with the MO, are found 
at g — r ~ 0.4, g ~ 19.5. The strong overdensity of stars at g — r ~ 1.5 are the local low mass dwarf stars. 

log[t/yrs]< 10.2), 30 % binary fracti on, a Salpeter Ini- 
tial Mass Function (Salpeter 1955), and three metal- 
licity bins, sufficient to describe the halo and th ick disc: 
[Fe/H1 = -0.7, [Fe/H]= -1.3, [Fe/H]= -2.2 [see |de Jong 
et al.| pilOt for 



O - g band 
x r band 



30 33 34 36 38 

Magnitude 

Fig. 3. — Example completeness data and fit for a field in the I 
= 150° stripe using the stars in the pointings which overlap. The 
completeness is determined by whether a star detected in either 
frame is recovered in the other frame. This is then binned by mag- 
nitude and presented as the percentage completeness of recovered 
divided by detected. Circles for the g-band objects and crosses for 
the r-band objects. 

stellar populations for comparison with t he data. 

In t he same manner as explained in |de Jong et al.| 
( |20T0| we used the S PSS g and r isochrones provided 



by Girardi et al. ( |2004 [). Given that both the thick disc 
and stellar halo are known to have old stellar popula- 
tions, we considered a fixed age range at ~13 Gyrs (10. 1< 



ate metallicities are inferred from the relative weight of 
these three template populations. The thin disc stars 
are avoided as they have a broader range of ages and 
metallicities which would make it difficult to disentan- 
gle from a combination of only three different metallicity 
templates. To avoid edge effects, the model templates 
were created for distance moduli between 7 (^250 pc) 
and 22 (~250 kpc) in steps of 0.2. 

The basis of MATCH is the Hess diagram, a binned 
CMD in which the value of each bin is the square root of 
the number of stars. Synthetic Hess diagrams are then 
created for a range of ages and metallicities initially as- 
suming 1 M /yr star formation rate (SFR) which is then 
scaled and combined to best match the observed dia- 
gram. The synthetic diagrams are convolved with the 
photometric errors and completeness profile of the data 
to provide a realistic comparison with the data. When 
comparing the observed CMD with the synthetic CMD, 
MATCH uses a Poisson Maximum Likelihood statistic to 
determine the best-fitting single model or linear combi- 
nation of models. 

We used stars in the magnitude range 18.5 < g < 
23.0 and 18.0 < r a < 23.0 and in the color range 0.1 
< 9o — f < 1-1 (Figure [4]). These color cuts ensure that 
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we do not include faint, red stars belonging to the thin 
disc, while the magnitude cuts secure that we do not 
include spurious objects such as misclassified galaxies. 
For each population template, MATCH provides the star 
formation rate in M Q /yr for each distance modulus bin. 
The star formation rates arc then converted into stellar 
mass density. 

Figure [4] shows an example of an observed CMD and its 
best-fit model CMD for a single frame in stripe I — 130°. 
For each single frame one Hess diagram is created. The 
observed CMD is shown on the left hand side; the region 
used in our analysis is depicted with the dot-dashed rect- 
angle. The four panels on the right hand side are: The 
observed Hess CMD (top left), the model Hess CMD 
(top right), the residual Hess CMD after subtracting the 
model from the data (bottom left) and the residual sig- 
nificance, based on the number of stars expected in the 
model Hess diagram (bottom right). The model CMD 
reproduces well the main features of the observed CMD, 
such as the main plume of old MSTO stars, especially 
since we assumed only a simple model population. 

3.2. Density and Metallicity gradients 

MATCH provides the star formation rate (SFR, in 
M Q /yr) corresponding to each population template for 
each distance modulus bin and this is transformed into a 
stellar mass density. Figure [5] shows the density profiles 
for the lower half of the I = 130° stripe plotted against 
heliocentric distance. A Sersic profile is then fit to the 
underlying stellar population and is shown overplotted 
on the data. The stellar density shows a clear decrease 
with distance and includes a conspicuous deviation in 
the distance range 7 < d (kpc)< 13. This "excess" in 
the density distribution, present in all our fields, is due 
to the MO. The steep exponential profile in the inner ~5 
kpc is due to the contribution of the thick disc popula- 
tion. From ~5 kpc onwards the inner halo component 
dominates up to ^20 kpc when the outer halo begins to 
rule, producing a flattening in the density profile. This 
correspond s well to the density profiles reported in de 
Jong et al |fl2010| ). 

The depth ot the data allows us to unequivocally delin- 
eate the density of stars in the MO. In order to obtain a 
clear detection of the MO and to find out if there are dif- 
ferences with height above the plane, we gathered all the 
fields corresponding to each stripe into two halves: upper 
and lower latitude, i.e., six halves in total, two per stripe. 
This improves the signal to noise of the MO as individual 
frames do not contain enough stars to perform the anal- 
ysis. To quantify the resultant overdensities, we removed 
the smooth background stellar density distribution and 
fitted a Sersic profile to the stellar mass density relation 
obtained from the CMD analysis. For each density pro- 
file, we took all the points within the 3a values of the 
Sersic fit and recalculated the density, thus removing the 
bulk Milky Way components from the distribution. The 
best-fit Sersic parameters can be found in Table [3] The 
resultant residual for each of the six fields is shown in Fig- 
ure [6j The upper panels represent the residuals of the 
totalmass density in the upper latitude set and the lower 
panels are the same for the lower latitude set. All of the 
residuals show a bump at the location of the MO. The 
number density is relatively constant across the stripes 
for each latitude range, however the lower latitudes are 



consistently denser than the higher latitudes. Figure |6] 
also shows the location of the MO, denoted by D, which 
was found by fitting a Gaussian profile to the residual 
density peaks. The line-of-sight depth of the MO rep- 
resents the full width at half maximum of the best fit 
Gaussian. The stellar number density of each of the six 
regions can also be found in Table [3j as well as the mean 
latitude, heliocentric distance and fine-of-sight depth. 

Given the similarities between each of the stripes, we 
gathered all the stripes together to obtain metallicity 
profile. The total mass-weighted mean metallicity profile 
is shown in Figure [7] Although the smooth underlying 
Milky Way population has not been subtracted, the MO 
is still clearly visible. The metallicity distribution is at 
slightly higher distances than seen in the density pro- 
files and deviates from the smooth background between 
9 and 14 kpc heliocentric. It reaches a peak metallicity of 
[Fc/H]~ —1.0 which is consistent with the photometallic- 
iti es of the MO as det ermined through SDSS photometry 
byllvezic et all (12008). 



4. DISCUSSION 

These deep CMDs of the MO at three different galactic 
longitudes (Z=130°, 150°, and 170°) and covering a range 
of galactic latitudes (+15° <6<+25°) allow us to accu- 
rately constrain the structural properties of the MO in 
these directions. Figures [2] [5] and [6]show that the MO is 
easily identifiable at all stages of the analysis: it appears 
as a strong main sequence type feature in the CMDs; it 
shows a clear excess above the Sersic fit to the bulk Milky 
Way components in the stellar density profiles; occupies 
a distinct distance range within the sensitivity limits of 
the method; and has a metallicity that strongly differs 
from the the background Milky Way population. Fig- 
ure [8] shows the locations of the MO with respect the 
Galactic centre and the Sun. A line has been drawn at 
the Galactic radius of 17.0 kpc for reference. Each de- 
tection is shown illustrating its dista nce unce rtai nty and 

we 



4.2 and 4.3 



the width of the feature. In Sections 4.1 
will discuss the various formation scenarios i n lig ht of the 
density profiles uncovered here. In Section |4.4[ we will 
discuss the implications of the metallicity finding and its 
relevance to the outer disc. 

4.1. The Monoceros Overdensity as a Tidal Stream 

We compare our results with the only two cu rrent 

pOOlal and 



numerical simulations of Martin et al. 



Penarrubia et al. (2005 1 . Figured shows the comparison 
between the two models in the <3alactic latitud e range 
from 5° to 2 5°, as probed by the survey. The |Martin 
et al. ( 2004a ) model (upper panels of Figure [9]) shows a 
slight decrease in Galactic latitude across the survey and 
describes a distinct stellar s tream predominantly below 
b = 20°. Penarrubia et al. (2005) (lower panels of Fig- 



ure [9]) has a tidal stream model which is found mostly 
at higher latitudes. In this manner, we should exp ect to 
see a decrease in density at higher latitudes for the |Mar-] 
tin et ah] (|2004a) model an d an increase in density for 
the [Penarrubia ct al. ( 2005[ ) model. In the middle upper 
and lower panels ot Figure |9l is seen that the observa- 
tions roughly match both models although the measured 
change in density (see TableJSl) is contrary to both mod- 
els. In terms of the height above the plane (right panels 



G 



CONN ET AL. 




1 2 0.5 I 




Fig. 4. — The left hand figure shows an observed CMD for a single frame of the stripe I = 130°. The long dashed rectan gle re presents the 
region used in our analysis. The four right hand figures show the output of the CMD fitting routine MATCH (see Section |3,1| , The upper 
panels show the observed CMD (left) and best-fit model to the observed CMD (right). The lower panels are residuals after subtracting 
the model from the observed CMD (left) and residual significance (right) based on a fit parameter including the number of stars in the 
observed and model CMD. The grey scale bar represents the number of stars in the Hess diagram bins. The analysis was made in the 
regions 18.5 < g a < 23.0, 18.0 < r < 23.0 and 0.1 < g — r a < 1.1 of the CMDs shown in Figure[2] In this way, we avoid the clump at 
g > 23.0 and g — r < 1.0 mainly due to unresolved background galaxies and local dwarf stars located at (g — r a ) > 1.5. As seen in 
Figure [3] we have 95% completeness at g = 23.0 and 90% at r = 23.0. 



TABLE 3 

Results from the MATCH analysis, Distance, depth and density. 



Region 


Mean latitude 
(6°) 


Distance Heliocentric 
(kpc) 


Depth of MO 
(kpc) 


Stellar Number Density 
(counts) 


Sersic Index / 
reduced \ 2 


Sersic Scale 
Length (kpc) 


130 Stripe (upper) 


24 


10.4±0.5 


3.8±0.4 


463 


7.0/1.4 


0.2 


130 Stripe (lower) 


18 


10.7±0.4 


3.7±0.4 


671 


8.0/1.2 


0.1 


150 Stripe (upper) 


21 


9.3±0.5 


1.8±0.2 


553 


13.0/1.2 


0.11 


150 Stripe (lower) 


17 


9.6±0.2 


3.1±0.4 


609 


6.0/1.5 


0.54 


170 Stripe (upper) 


25 


10.2±0.3 


2.5±0.3 


427 


5.15/1.5 


2.43 


170 Stripe (lower) 


20 


10.5±0.2 


1.8±0.2 


588 


3.55/1.6 


0.47 



of Figure J9| the overdensity seems to be slightly closer 
at higher latitudes than at lower latitudes, although with 
the errors it is consistent with a vertical feature. Both 
models seem to bracket a possible tidal stream scenario 
for the MO as determined through this survey. Although 
the structure of the stream seen in the data is not com- 
patible with either simulation it is difficult to exclude a 
tidal stream solution since the large number of param- 
eters practically ensure a suitable model is likely to be 
found. 

4.2. The Monoceros Overdensity as the Galactic Flare 

The MO is a low-latitude stellar structure and, as such, 
could be related to the generic structure of the disc. Al- 
tho ugh many investigations have pursued this possibil- 
ity (iMomany et al 112004} jMoitinho . et al l|2006l |Momany 
|et al.|2l)U6| |L6pez-0orredoira et al.|2007 frTammersley fe 



Lopez-Corredoira 2011 1 the distance to the MO typically 
precludes a definitive conclusion since the MO stars are 
faint and removing contaminants is highly proble matic. 



Recently, |Hammersley & Lopez-Corredoira |(|2011|), have 



attempted to show that the stellar profiles seen in the 
directions of the MO are compatible with the flaring of 
the galactic disc. The flare is described such that beyond 
a certain radius, the disc rapidly thickens and becomes 
prominent above the plane, replicating the effect of the 
MO stars. They sample a small range of galactic longi- 
tudes, mostly in regions unaffected by the galactic warp, 
and fit a small range of flare models to the SDSS CMDs. 
They conclude that the stellar counts can be accounted 
for with the galactic flare starting at 16 kpc galactocen- 
tric and using a scale length of ~4.5±1.5 kpc. 

Although the CMD fitting method presented here dif- 
fers significantly from the star count approach used in 
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Fig. 5. — Mean stellar mass density distribution for the stripe 
with I = 130°, b < 21°. The stellar density decreases with increas- 
ing distance with a conspicuous deviation at around 10 kpc from 
the Sun, corresponding to the location of the MO. A Sersic profile 
is fitted to the density distribution to remove the Milky Way back- 
ground and allow for a cleaner detection of the MO. Error bars 
were calculated from the uncertainty associated to the density of 
each bin from the Monte Carlo tests. This is achieved by mul- 
tiply resampling the stars within the CMD space and re-running 
MATCH on each resampled population. 



Fig. 6. — The residual density of the stream after subtracting a 
Sersic profile fit, in the distance range of 1 — 25 kpc, to the smooth 
Milky Way component. The top panels are the residuals from the 
higher latitude half of each stripe, while the bottom panels are the 
residuals from the lower latitude half of the stripes. See Table|3]for 
the obtained values from this analysis. The distance to eacnMO 
detection is determined through a Gaussian fit to the residuals. 
The width of the MO is the FWHM of that Gaussian fit. This is 
shown as 'D' and 'AD' in each sub panel. The error on the distance 
is related to the error on fitting the Gaussian to the residual density 
profile. 



TABLE 4 

Flar e parameters used in Equations [7] 
from ffammersley & lopez-corredoira 
( |2011| . 



Parameter 



Value 



thin disc scale height {h z ^hin,o) 186 pc 

thin disc scale length (h,R thin) 2400 pc 

thick disc scale height {h z thick, o) 631 pc 

thick disc scale length (Iir thick) 3500 pc 

Solar radius (i? ) 7900 pc 



Hammersley & Lopez-Corredoira (2011), we have fitted 
their hare models to our dataset across a large parame- 
ter space of flare scale lengths and flare onset positions 
to further investigate this scenario. In this regard, we 
also use the equations below to define the flare as de- 
scribed in their paper. Table [4] lists the constants used 
in the model. For convenience, we reproduce them here 
(Equation [2]) , note that Ri is the radius at which the 
flare starts; A is a scale factor for the density; pthin is 
the density of the thin disc; pthick is the density of the 
thick disc; phaio is the density of the halo; h r f is the flare 
scale length; R is the galactocentric radius and z is the 
height above the disc. 



= A 



Ptotal — Pthin + Pthick " 
R - R @ 



Phalo 



(2) 



fl»,tto(fl) 



exp 



exp 



^z,thin(fi) 



Pthick = 0.09A 



hzA 



exp 



R - Rr. 



exp 



Phalo = 1.4 x 1(T 3 A 



exp [10.093 (l-(^) )] 



^z, thin/thick 



(R) 



(R sp /R Q y/s 

^-z,thin/thick,0j < Ri 



2,thin/thick 



R sp = VR 2 + 2.52z 2 



, exp(^),i?>R i 



To test this model against our MO density profile we 
have varied both the onset point of the flare and its scale 
length. Figures 10 and 11 show the resulting chi-square 
space with 1, 2 and 3cr contours for the global and in- 
dividual fits to the data, respectively. The flare onset 
position varies from 1 to 21 kpc and the scale length has 
been varied from to 10 kpc, both of which were iter- 
ated in 100 pc steps. Each model was compared against 
the data and the chi-square value was determined for 
each point in the parameter space. Flare models with 
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Fig. 7. — Mass-weighted mean metallicity from the MATCH fits, 
averaged over all latitudes and longitudes. At distances less than 
10 kpc, the metallicity gradient is steep due to the thick disc halo 
transition. The MO, around 10 kpc, also appears to be a distinct 
feature in the mean metallicity compared to the 'underlying' dis- 
tribution seen again beyond ~15 kpc. Below 7 kpc, the stellar 
populations are most likely too complex for our simple approach 
as seen in the large error bars irregular profile. F or a complete 
discussion on the metallicity of the MO, see Section |4.4| 



TABLE 5 

Best Flare Model Parameters determined by the 
minima of chi-square map as seen in Figure |TT| with 

THE UNCERTAINTIES DETERMINED BY THE la CONTOUR 
LINE. 



Field Onset Position (kpc) 


Scale Height (kpc) 


Global 


12 6 +0 - 3 


9 i +0.3 

z,A -o.a 


130 (Upper) 


11.7±?» 


3 2+ 1 - 4 


130 (LowerF] 


12.8 


1.7 


150 (Upper) 


13.2+0J 


1 5+ - 7 


150 (Lower) 


13 5+ - 7 


1 8+ ' 6 


170 (Upper) 


11 3+ 11 


3 5+ 1 - 
°*°-1.3 


170 (Lower) 


13 -°±i:i 


9 9+I.O 
z - z -0.7 



a The lower field in the 130 stripe has no error estimate 
since the min ima used here is not the absolute minima found 
in Figure [TT| 



small scale lengths describe sharp features in the density 
profile while large scale lengths have long slowly varying 
density profiles. The data has been fitted both to find a 
global solution for the flare considering all the data and 
indivi dual ly to highlight the differences between fields. 
Figure [Tl] shows that the MO density profile typically re- 



FlG. 8. — Galactocentric locations and distances of the MO in 
each of the three stripes. The thick line shows the error on the 
distance estimate while the thin line shows the depth of the stream. 
The dashed line is the direction of each of the stripes. The lower 
latitude fields plotted to the right of the line and the higher latitude 
fields to the left. The Sun is located 8 kpc from the centre and the 
solid line is a galactocentric circle with radius 17.0 kpc. The lines 
from the galactic centre are visual aids to highlight the circularity 
of the curved line. 
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Fig. 9. 



L eft Panels: [Martin ct al. (2004aJ (top) and |Pefiarrubia| 
|et al.| (["2005} (bottom) m tne same region as surveyed m tms pa- 
per. The green diamonds show the mean latitude of the upper half 
of this survey and the blue diamonds show the mean latitude of 
the lower half of this survey. Centre Panels: show the top down 
galactocentric view of the models in the Second Quadrant. The dis- 
tances to MO are shown with both green and blue filled diamonds. 
Right Panels: show the cross-section view of the tidal models with 
galactocentric radius on the x-axis and height above the plane on 
the y-axis. 
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Fig. 10. — Chi-square space after fitting for a global solution to 
the galactic flare model. The onset position and scale length arc 
free parameters in the fit with the contours showing the 1, 2 and 
3cr deviations. 



quires the flare model to have a very short scale length, 
with a global solution of 2.1 kpc and an onset radius 
of 12.6 kpc (see Table [5] for the best fit in each field). 
In Figure [12] the globalmodel has been overplotted on 
each of the MO density profiles for comparison with the 
non-flare model shown as a solid-blue line. The model of 
the flare used here is very basic and so does not include 
a prescription for known Galactic features such as the 
Warp. This is clearly evident in the fields closer to I — 
90° where the Warp becomes stronger. The mismatch at 
small heliocentric distances for the lower half of I = 130° 
stripe is an example of this. 

In general, this basic flare model can be fit to the data 
within the uncertainties for the majority of the points. 
While the presence of the warp is clearly responsible 
for the discrepancies at small heliocentric distances, it 
is unclear whether the differences between the MO den- 
sity profile and the generic properties of the flare model 
should be explained by the noise in the data or an intrin- 
sic irregularity of the outer disc. The flare parameters 
fou nd he re are consistent with that found by |Mateu et| 
al. (2011 ) (onset ~ 11.5 kpc, scale length ~1.6 kpc), using 
KR Lyrae stars to trace the outer thick disc. However, 
this is a much shorter scale length and onset radius than 
that found by Hammersley & Lopez-Corredoira (20111 
(onset ~ 16 kpc, scale length ~4.5±1.5 k pc). This differ- 
ence potentially arises fro m the fact that [Hammer sley fc| 
Lopez-Corredoira (2011) push to stars with increasing 
photometric e rror s which will inherently smear out their 
result. Figure [TT| shows that a scale length of ~4.5±1.5 
kpc is feasible at the 2 — 3cr level but a flare onset position 
of 16 kpc is not likely with our data. Improved number 
statistics might be necessary to reduce the uncertainties 
in the MO density profile and thus we could determine 
whether the deviations from the smooth model can be 
considered significant or simply the nature of the outer 
disc itself. 



There is the possibility that the MO could be explained 
through a disrupted disc scenario whereby the disc inter- 
acts with a massive dark matter sub-halo. In this sce- 
nario, no new stars are added to the disc but rather the 
existing disc stars are swept or migrated into large spiral 
or ring-like stru ctures. Models illustrati n g this scenario 
can be found in Kazantzidis ct al. (2008); Younger et al. 
( p008l l ; |Purcell et al.| ( |201l| ); |G6mez eFaL]fl201l| i. All the 
authors find that the structures are typically ~4 Gyrs old 
and so are relatively long-lived. 

Distinguishing between these and a tidal stream is 
difficult and most likely requires detailed velocities or 
chemical abundance information. Fortunately, some in- 
sights into whether this scenario is feasibl e can be seen 
in the stellar density profiles as shown in |Purcell et al 
(2011 1. In their model, the resulting stellar density pro- 
file with heliocentric radius has significant substructure. 
Since the overdensity is not created with new stars but 
rather a rearrangement of the disc, it follows that cre- 
ating an overdensity naturally produces a corresponding 
underdensity. In this manner, both their light and heavy 
Sagittarius-like dwarf galaxy encounters induce signifi- 
cant underdensities in the disc, adjacent to the ring-like 
overdensity o f arou nd 1.2 dex (see Fig 4b and S7 from 



Purcell et al. (2011)). Crucially, there is no evidence for 
these underdensities in our dataset that could possibly 
match the dramatic change in the stellar density profile 
as suggested by their model. Additionally, to place stars 
at the location of the MO detections ~5 kpc above the 
plane, is limited to their simulation with a heavy Sagit- 
tarius dwarf galaxy. The light version is unable to have 
such an impact on the scale heights of the disc stars. 



4.4. The Metallicity of the Monoceros Overdensity and 
the Outer Disc 

The metallicity of the MO has been measured both 
photometrically and spectroscopically with a variety of 
results. Photometrically, the metal licity has bee n found 
to be [Fe/H] 0.95 ± 0.15 from |Iyezic et al l ( |2008[ ), 

[Fe/H] ' '" " 



d201l|) and \WJEp 
t has been reportc 



-1.5 ± 0.1 from pesaretal 
T.0, this study. Spectros copicallyTit has been reported 
as [Fe/H]~ -1.6 ± 0.3 by JYanny et al] ( [20031 ), [Fe/H]~ 

- 0.4±0.3 by|Crane et al.| ([2003|),'<|Fe/Hj> 1.37 ±0.04 

by|Wilhelm et al.| ( |2005[ ) =TTM < [ Fe/H]< -0.1 by|Chou| 
|et al.| ( |2010|) and most recently by |Meisner et al.| ( f27JPZ| ) 
with [Fe/H] = —1.0. From this, it is clear that the 
MO is a complex stellar population that is consistently 
metal- poor. 

Characterizing the disc at the distances of the MO is 
difficult and so the expected metallicity profile needs to 
be extrapolated from our understanding at smaller galac- 



outer disc 1 


Co§kunoglu et al. 


2012| Lee et al. 2011 


Cheng 


et al. 2012 Ivezic et al. 2008 


Vong et al.|2005| to 


under- 



SEGUE survey, Lee et al. (2011) traced the metallicity 
of the thin and thick discs beyond 2 kpc from the Sun 
showing that the mean [Fe/H] o f the thin dis c is ~ —0.2 



4.3. The Monoceros Overdensity as a Perturbed disc 



and for the thick disc ~ —0.6. Cheng et al. (2012) also 
finds similar results with SEGUE finding that stars be- 
tween 1.0 - 1.5 kpc above the plane have a metallicity of 
0.3 <[Fe/H]<-1.0, centered on [Fe/HK-0.6 In terms of 



metallicity gradients, Coskunoglu et al. (2012) has shown 
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130 Stripe (Upper) 150 Stripe (Upper) 170 Stripe (Upper) 
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Fig. 11. — Each panel shows the chi-square space after fitting a galactic flare model (Eqn.[2l for the upper and lower halves of the three 
stripes of the survey with the global solution shown in the top-left panel. The figures show tne range of onset points for the galactic flare 
and its the scale length for the model tested with the resultant chi-square plotted in color. The contours shown in each sub-panel delineate 
the 1, 2 and 3a lines. In general, the best-fit models are exclusively with very small scale lengths and the onset radii increases from ~12 
kpc at / = 130° to ~14 kpc at I = 170° inline with the 2 kpc change in distance seen in CMD fitting analysis. The small scale length, in 
particular, reveals the nature of the MO to be one of a short transient feature and not an extended, generic component of the disc. The 
minima at onset values less than 10 kpc and large scale lengths seen in the lower halves of the I = 130° and 150° stripes are due to the 
model not having a prescription for the warp and the densities mismatching at small heliocentric distances. 



with the RAVEE3 dwarf stars that while the thin disc 
decreases in metallicity by —0.043 dex/kpc -1 , the thick 



disc i s essen tially fiat. Both Lee et al 
et al. ( 2012 ) also find the thic k disc to 
gradient 



(|201_Lj) and [Cheng 
rave no metallicity 



Yong et al. (20051 explored the outer disc in 
the third galactic quadrant using open clu sters and also 
found a flat distribution with [Fe/H]~-0.6. Bensby et al 



(2011 ) confirm the metallicity gradient in the thin disc as 
their target stars in the galactocentric distance range of 9 
- 13 kpc have thin disc abundance patterns with a mean 
metallicity of [Fe/H]~ — 0.48±0.12, which is significantly 
more metal-poor than the local thin disc stars. Extrapo- 
lating the metallici ty gradient of —0.066 dcx/kpc" 1 from 
Cheng et al. ( 2012 ) to the radii of the MO, we find a pre- 
dicted thm disc metallicity of [Fe/H]~ —0.44, which is 
still more metal rich than all estimates of the MO metal- 
licity. 

A final possibility remains that the metallicity derived 

10 RAdial Velocity Experiment 



through isochronc fitting is wrong, simply because we 
have utilized old metal poor isochrones which would be 
unsuitable for a thin disc population. If the MO were 
thin disc stars then a 4 Gyr isochrone would better repre- 
sent such a population. In this case, the isochrone would 
be bluer by (g — r a ) ~ 0.1 and so could feasibly be con- 
sistent with our data. The difficulty with this approach is 
that a 4 Gyr old main sequence turn-off star is at least 1 
magnitude brighter than the corresponding 10 Gyr star. 
At the turn-off magnitudes seen in the data, g ~ 19, 
this translates into an additional 5 kpc in line-of-sight 
distance, placing the MO at ^21.5 kpc Galactocentric. 
Naively extrapolating the metallicity of the disc to these 
distances results in [Fe/H]~ —0.7. The MO though, is 
now 10 kpc above the Plane and so at each turn it be- 
comes harder to associate thin disc stars with the MO. If 
the thick disc truly exhibits no change in its metallicity 
distribution with radius then the MO also remains more 
metal poor than the thick disc at these large radii. 
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Fig. 12. — CMD fitting analysis density profiles for each of the upper and lower halves of the three stripes of the survey overplottcd 
with the global galactic flare model (green points) plotted against heliocentric distance. The non-flare model is shown as a solid-blue line. 
A discrepancy within the first 5 kpc is due to the presence of the warp which is unaccounted for by this model and the I = 130° lower 
half-stripe has a minima due to this mismatch. Overall, the global flare model is within the error bars although in several fields it does not 
trace the peak and is consistently lower than the data at large galactic radii. 



Given this understanding of the outer disc, we can in- 
terpret the likelihood of the different formation scenarios 
with the metallicity finding of this study. It is important 
to note that this method of determining the metallicity 
relies on the bulk properties of the stars in the CMD and 
is not a direct measure of distinct components like the 
Thick disc and Halo. Rather, the profile as seen in Fig- 
ure [3 shows how the contributions of the Halo stars be- 
come more dominant with increasing distance and so the 
average metallicity of the stellar populations present is 
increasingly more metal poor with heliocentric distance. 
The MO is therefore a distinct population which abruptly 
appears against this smooth transition to a pure Halo 
population beyond the disc. 

Tidal stream scenario: Distinguishing between local 
disc stars and stream stars from a merger is perhaps 
clearest in the ch emical abundance pa tterns as shown 



in the review by Tolstoy et al. (20091. There are dis 



tinct chemical differences between local MW stars and 
stars from nearby galaxies that reveal their different en- 
richment histories. Recent studies of the M O using spec- 
troscopically determined abundances (e.g. Chou et al. 



( |2010| ; |Meisner et al.| ( |2012[ )), show that the chemical 
properties of the MO are closer to a Large Magellanic 
Cloud or Sagittarius Dwarf galaxy type abundance pat- 
tern than a pure Milky Way disc population. This offset 
in metallicity between the MO and outer d isc suggests 
poten t ially a different origin for these stars. Penarrubia 



et al.| ( 2006 ) suggest the outer disc could have been ere 
ated through a series of mergers in which case the abun- 
dance pattern and the consistently metal-poor nature of 
the MO me mber stars are supp ortive of this scenario. 
Additionally, Carollo et al. (2010) discuss the similarities 
between the metal- weak thicFdlsc (MWTD) and the MO 
suggesting the two may be related. Indee d the MWTD 
itself is presented in Carollo et al.| ( 2010[ ) as distinct to 
the canonical thick disc and as such is possibly the result 
of a merger with the Milky Way disc. Together the evi- 
dence builds that the MO is an accretion event although 
there is no viable progenitor and its passage through the 
outer disc is still unknown. 

Galactic Flare scenario: Although the flare model of 



Hammersley & Lopez-Corredoira (2011) does not make 
any predictions about the metallicity ot the stars, our un- 
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derstanding of the disc can be used to determine whether 
the MO metallicity is consistent with a Milky Way pop- 
ulation. It is clear from Figure[7]that the stars along our 
lines of sight have a steadily declining metallicity with 
distance and the stars bracketing the MO typically have 
abundances of -1.2 < [Fc/H]< -1.5. Thus the MO ap- 
pears distinct in the outer disc as more metal-rich than 
the nearby stars. A compari son between these metallic- 
ities and those described in Ivezic et al. (2008) suggest 
that at these distances we are beginning to probe the in- 
ner halo prior to MO and beyond the MO there is a clean 
halo sample. Clearly, these stars are apart from the main 
disc population but it is difficult to explain why the MO 
is so metal poor if it is simply an extension of the under- 
lying disc. The flare is undoubtedly a real phenomenon 
but to what extent and what influence it has in the outer 
disc is uncertain. 

Perturbed disc scenario: The star s which are per t urbed 
into the MO-like structure seen in Purcell et al. (20111 
are sourced from across the entire disc. The member 
stars are migrated from inside and outside of the final 
location and so the resultant metallicity should be an 
average of these contributing locations in the disc. Since 
the disc, in general, is more metal rich than the MO and 
there are very few locations within the disc which could 
supply stars more metal poor than the MO, it is highly 
unlikely that an aggregate population as proposed by this 
model, could achieve the metal poor status of current 
set of MO metallicity estimates. Since our findings too, 
confirm the metal poor nature of the MO, the perturbed 
disc scenario with its pred i cted o bservable properties, as 
described by Purcell et al. (2011 1 is not feasible given the 
data. 



5. CONCLUSION 

We have presented new distance, density and metallic- 
ity measurements for the stellar Monoceros Overdensity 
(MO) in the outer Milky Way, based on SUPRIME-CAM 
wide field imaging data and a CMD fitting analysis. Our 
distance measures are the most quantitative estimates to 
date for the MO. 

The MO appears as a wall of stellar material at 
roughly 10 kpc from the Sun at the galactic longitudes 
of 130°, 150° and 170°, and galactic latitudes of +15° < 
b< +25°. Detections of the MO have been confirmed 
between 3-5 kpc above the plane and consist of a 
metal-poor population with an average metallicity of 
[Fe/H]~ -1.0. 

We consider these findings in the light of the three for- 
mation scenarios currently in the literature: (i) a tidal 
stream origin; (ii) the galactic flare; and (iii) the per- 
turbed disc. We find that: 

(i) Tidal stream models from the literature bracket the 
distances and densities we derive for the MO. Fur- 
thermore, recent results for the chemistry of stars 
in the MO support an extragalactic origin. This 
suggests that a tidal stream model can be found 
that would fully fit the MO data. On the other 
hand, the large parameter space available for this 
model: the orbit, mass, inclination and eccentricity 
of the merger, amongst others, presents the danger 
that such a fit - while possible - might not be the 
true explanation for the MO. 



(ii) The flaring of the galactic disc provides another 
possibility for explaining the presence of these stars 
at large distances from the plane. We fitted a 
large range of galactic flare models finding a so- 
lution with a mean onset radius of 12.6 kpc and a 
scale length of 2.1 kpc that is a reasonable match 
the data. Th is is similar to the findings of |Mateu| 



et al. (20111 bu t is much smaller than the mod- 
els sug gested by Hammersley fc Lopez-Corredoira 
(2011). The main difficulties with the flare model 
are: (a) whether the basic flare model used here 
while consistent with the data would be applica- 
ble across wider latitude and longitude ranges and 
(b) the metallicity ([Fe/H]) derived in this paper as 
well as the determinations from other sources (see 
Section 4.4) are building a consistent picture that 
the disc is too metal-rich to source the MO stars. 
If the disc can be shown to be metal-poor at these 
radii then the flare scenario is indeed a possibility. 



(iii) The perturbed disc scenario makes clear testable 
predictions about the metallicity and stellar den- 
sity profile of a MO-like feature. Both of these are 
incompatible with the data: the MO stellar density 
profile does not contain the significant underdensi- 
ties predicted by the model while the metallicity of 
the MO is too metal-poor even for a population of 
stars sourced from across the disc. 

It is clear that the MO still lacks the observational ev- 
idence required to unequivocally determine its origins. 
However, the deep observations we have presented here, 
coupled with CMD-fitting techniques are able to con- 
strain its properties to much greater precision than has 
previously been possible. We have ruled out the 'per- 
turbed disc scenario' for the MO, and found key prob- 
lems that must be solved if the MO is to be explained 
by a flared disc. Given the distance to the MO and the 
uncertainty over its origin, it is crucial to minimize the 
photometric errors so as to limit their impact when de- 
riving its properties. Further studies of the MO should 
include high precision photometry to better constrain the 
physical dimensions of the MO coupled with high resolu- 
tion spectroscopy for a detailed abundance analysis. This 
combined approach seems best suited to unravelling the 
origin of the MO feature. 
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TABLE 6 
Column information and format. 



i Number 


Label 


Type 


Column Number 


Label 


Type 


1 


Right Ascension (Hours) 


Integer 


11 


Y pixel g 


Real 


2 


Right Ascension (Minutes) 


Integer 


12 


g mag 


Real 


3 


Right Ascension (Seconds) 


Real 


13 


g mag error 


Real 


1 


Declination (Degrees) 


Integer 


11 


g classification 


Integer 


5 


Declination (Minutes) 


Integer 


15 


X pixel r 


Real 


6 


Declination (Seconds) 


Real 


16 


Y pixel r 


Real 


7 


Galactic Longitude (I) 


Real 


17 


r mag 


Real 


8 


Galactic Latitude (b) 


Real 


18 


r mag error 


Real 


9 


Chip 


Integer 


19 


r classification 


Integer 


10 


X pixel (g) 


Real 


20 


E(B-V) 


Real 



CAT (http://www.starlink.ac.uk/topcat/) was used ex- 



tensively in the preparation of this paper. 
Facilities: SUBARU (SUPRIME-CAM) 

APPENDIX 
ONLINE DATA 

The data presented in this paper is being made available online and all issues related to the data can be addressed to 
conn@mpia-hd.mpg.de. The data consists of all objects extracted from each reduced frame in the survey, galaxies and 
stars, although poor data has been excised from the final catalogue. Individual reduced frames will be available upon 
request. The catalogue mostly follows the layout of the Cambridge Astronomical Survey Unit Wide Field Camera 
Pipeline (Irwin & Lewis |2001 ) with the addition of the Galactic coordinates for each object, Right Ascension and 
Declination are in J 2000.(1 

The descriptors for each column in the catalogue is shown in Tablel6]and an example of the data is shown in Table [7J 
The classification scheme is as follows: Stellar are —1, Possible Stellar are —2, non-Stellar/Galaxy are 1, Noise is 0, 
Possible non-Stellar/Galaxy is —3, Crossmatch problem is —8, Saturated object is —9. The entire catalogue contains 
3.4 millions objects and the breakdown per filter in the entire catalogue is —643,000 Stellar, —550,000 Possible Stellar, 
-747,000 non-Stellar/Galaxy, -411,000 Possible non-Stellar/Galaxy in the 5-band. The r-band has -625,000 Stellar, 
-432,000 Possible Stellar, -764,000 non-Stellar/Galaxy, -320,000 Possible non-Stellar/Galaxy. 

The data here has been extinction correcte d following the s tanda rd prescription when t he ex tinction is less than 
E(B — V) < 0.1, otherwise the relation from Bonifa cio et al. ( 2000[), as shown in Equaton |Al| has been used. The 
dust values have been extracted from the |Schlegel, Einkbeiner, & Davis ( |1998 ) maps using the dust_getval.c program. 



dust value = 0.1 + 0.65(E(B - V) - 0.1) (Al) 

The data has also been corrected for airmass. The airmasses and field centers for each field can be found in Table [51 
Individual objects do not have an identifier which relate them to a particular field however most objects will be easily 
matched with its corresponding field center. Objects in overlap regions can be associated with a particular field based 
on the chip in which they reside. In this regard, users should note that duplicate observations of the same object have 
not been removed from the catalogue. They have been left in to allow the user to gauge the relative depths of each 
pointing. 
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TABLE 7 

Example data from the catalogue 



v_'01 1. 


rv.i o 
K^Ol A. 


OOl o. 






OOl 0. 




L^Ol o. 


ooi y. 




5 


32 


29.55 


62 


16 


34.7 


149.930191 


15.2970066 


i 


1553.4 


5 


32 


40.61 


62 


14 


34.4 


149.971298 


15.2990808 


i 


1162.8 


5 


32 


30.45 


62 


7 


57.5 


150.060394 


15.2272711 


i 


1517.3 


5 


32 


30.43 


62 


12 


27.5 


149.992889 


15.2644424 


i 


1522.6 


5 


32 


42.7 


62 


16 


15.0 


149.948212 


15.3164396 


i 


1087.9 


5 


33 


7.83 


62 


6 


31.6 


150.119278 


15.2787237 


i 


184.5 


5 


32 


35.22 


62 


5 


13.7 


150.106125 


15.2127686 


i 


1343.1 


5 


32 


58.89 


62 


14 


42.0 


149.987595 


15.3309498 


i 


512.7 


5 


33 


6.82 


62 


6 


24.8 


150.119965 


15.2760792 


i 


220.4 


5 


32 


41.21 


62 


15 


37.5 


149.956116 


15.3087721 


i 


1141.0 


Col 11. 


Col 12. 


Col 13. 


Col 14. 


Col 15. 


Col 16. 


Col 17. 


Col 18. 


Col 19. 


Col20. 



3848.9 


22.552 


0.018 


-1 


1550.4 


3851.3 


21.162 


0.011 


-1 


0.113 


3248.3 


22.56 


0.019 


-1 


1159.7 


3250.6 


21.157 


0.011 


-1 


0.145 


1251.9 


22.555 


0.019 


1 


1514.4 


1254.1 


20.76 


0.0080 


1 


0.126 


2602.9 


22.527 


0.019 


1 


1519.6 


2605.2 


21.508 


0.015 


1 


0.102 


3758.9 


23.611 


0.019 





1084.8 


3761.4 


3.616 


0.0 





0.198 


840.7 


22.586 


0.019 


1 


181.2 


842.9 


20.898 


0.0090 


1 


0.122 


436.8 


22.496 


0.019 


-1 


1340.2 


438.9 


22.293 


0.03 


-1 


0.114 


3300.5 


22.61 


0.019 


-3 


509.4 


3302.9 


21.231 


0.011 


1 


0.109 


806.1 


22.596 


0.019 


-1 


217.2 


808.2 


21.087 


0.01 


-1 


0.108 


3567.9 


22.613 


0.019 


-1 


1138.0 


3570.4 


21.07 


0.01 


-1 


0.109 
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TABLE 8 

Field centers for the individual pointings of the survey, g— and t-band observations have the same pointing centers. 



Field Name 


Galactic Longitude (I) 


Galactic Latitude (6) 


R..A. (degrees) 




Dec (degrees) 


Airmass (g] 


) Airmass 


130. 


.01 


129. 


.505203 


15 


.2144403 


41.34738719353522 


76 


.60936674102940 


1 


.83 


1 


.88 


130. 


02 


129. 


.498260 


15 


.5623217 


41.99728376377872 


76 


.92390791959453 


1 


.85 


1 


.89 


130. 


.03 


129. 


520309 


15 


.9293108 


42.83053801772657 


77 


.24082916240711 


1 


.87 


1 


.91 


130. 


.04 


129. 


.490295 


16 


.3293571 


43.57344637209229 


77 


.60785831085641 


1 


.88 


1 


.93 


130. 


.05 


129. 


.477982 


16 


.7765141 


44.54271403857756 


78 


.00559018571801 


1 


.90 


1 


.95 


130. 


.06 


129. 


488922 


17 


.2290936 


45.68634624732895 


78. 


.39321646040258 


1 


.93 


1 


.98 


130. 


.07 


129. 


.481201 


17 


.7492542 


47.01153201618068 


78 


.84294581324289 


1 


.95 


2 


.00 


ISO- 


08 


129. 


.486389 


18 


.2244434 


48.36961137572538 


79 


.24196293987690 


1 


.98 


2 


.03 


ISO. 


.10 


129. 


.891327 


18 


.6489811 


51.37118049084526 


79 


.37541996407468 


1 


.99 


2. 


.05 


130. 


.11 


130. 


.495422 


18 


.4339542 


53.12960548181101 


78 


.86147318029894 


1 


.96 


2. 


.03 


130. 


.12 


129. 


.932236 


19 


.6258717 


54.88780147544033 


80 


12706264257046 


2 


.04 


2 


.10 


130. 


13 


130. 


.272293 


19 


.6864243 


56.53253062677564 


79 


.96761144072845 


2 


.03 


2. 


.10 


130. 


.14 


130. 


.025726 


20 


.5192966 


58.82546202768954 


80 


.74218057326074 


2 


.08 


2. 


.16 


130. 


.15 


130. 


352402 


20 


.5899010 


60.48630003355603 


80 


,58115326477444 


2 


.08 


2. 


.15 


ISO- 


16 


129. 


.884949 


21 


.1346321 


60.96542286637656 


81 


.27529378945327 


2 


.12 


2. 


.20 


ISO. 


.20 


129. 


500458 


23 


.1141567 


70.27297710747496 


82 


.82702811918269 


2 


.20 


2 


.21 


130. 


.21 


129. 


.786774 


23 


.3386803 


72.92685925574828 


82 


.73512051347623 


2 


.19 


2 


.21 


130- 


-22 


129. 


.711731 


23 


.7662678 


75.57591317062941 


83 


.01743603239895 


2 


.21 


2. 


.23 


130. 


.23 


129. 


709274 


24 


.1529255 


78.39220828275825 


83 


.20574701910283 


2 


.23 


2. 


.26 


130. 


.24 


130. 


.074417 


2d 


.3651276 


81.19850907262055 


83 


.00000442870187 


2 


.21 


2 


.25 


130. 


.25 


130. 


.503357 


24 


.6074963 


84.23756455977386 


82 


.73757597853790 


2 


.20 


2. 


.24 


ISO- 


26 


129. 


883682 


25 


.2132816 


87.35185357260404 


83 


.47482960191130 


2 


.26 


2. 


.30 


ISO. 


.27 


130. 


.162216 


25 


.6129627 


91.31864264785848 


83 


.34836976965003 


2 


.25 


2 


.31 


130. 


.28 


130. 


.274887 


25 


.7499065 


92.67235820452719 


83 


.28326763729932 


2 


.25 


2 


.31 


130. 


29 


129. 


845047 


26 


.3427334 


97.17409024268696 


83 


.77700149427237 


2 


.30 


2 


.37 


150. 


.01 


149. 


.367172 


14 


.7309074 


81.47717324443296 


62 


.45923524113733 


1 


.54 


1 


.75 


150. 


.02 


150. 


.021484 


15 


.0184927 


82.70019456402970 


62 


.06101665160047 


1 


.53 


1 


.44 


150. 


03 


150. 


568130 


15 


.1846695 


83.56255315448394 


61 


.68566109101802 


1 


.53 


1 


.43 


150. 


04 


149. 


.267120 


15 


.8481741 


83.48413227595164 


63 


.10354272409283 


1 


.55 


1 


.46 


150. 


.05 


149. 


.467026 


16 


.1752510 


84.32222788603623 


63 


.09321359738084 


1 


.55 


1 


.46 


150. 


.06 


149. 


739899 


16 


.3192616 


84.87544284777340 


62 


.93090358029885 


1 


.54 


1 


.46 


150. 


-07 


149. 


.869659 


16 


.5827961 


85.51538476758475 


62 


.94394993994977 


1 


.54 


1 


.46 


150. 


.08 


150. 


.418961 


16 


.8561764 


86.57224144450691 


62 


.60087744573428 


1 


.53 


1 


.45 


150. 


.09 


149. 


855774 


17 


.7393131 


87.79741617372791 


63 


.47297203588319 


1 


.55 


1 


.47 


150. 


10 


149. 


978928 


17 


.9487495 


88.33340417316249 


63 


.45655395917793 


1 


.55 


1 


.47 


150. 


.11 


149. 


.595703 


18 


.8632870 


89.88188730710624 


64 


.16482959621780 


1 


.58 


1 


.50 


150. 


.12 


149. 


404037 


18 


.9116516 


89.81690579089114 


64 


.35037968901068 


1 


.57 


1 


.50 


150. 


13 


149. 


.469101 


19 


.4710522 


91.06547830802131 


64 


.51266042943159 


1 


.58 


1 


.50 


150. 


-14 


149. 


.498032 


19 


.8628922 


91.93500344111298 


64 


.63475959071147 


1 


.58 


1 


.51 


150. 


.15 


149. 


465454 


20 


.0847473 


92.39163688508779 


64 


.74451585942646 


1 


.58 


1 


.51 


150. 


.16 


149. 


.472748 


20 


.4829426 


93.27157679386292 


64 


.88004096575455 


1 


.59 


1 


.52 


150. 


-17 


150. 


229660 


20 


.8936558 


94.73593492322217 


64 


.35456695208589 


1 


.58 


1 


.51 


150. 


.18 


149. 


.338928 


21 


.2984161 


94.99382608269545 


65 


.27239075985462 


1 


.60 


1 


.53 


150. 


.19 


150. 


.334595 


21 


.2837868 


95.66243700857011 


64 


.38920170930041 


1 


.58 


1 


.51 


150. 


-20 


150. 


471024 


21 


.6856174 


96.63736873754850 


64 


.39421276382994 


1 


.58 


1 


.51 


150. 


.21 


149. 


.105026 


22 


.5320530 


97.67760865552893 


65 


.85253465676418 


1 


.62 


1 


.55 


150. 


.22 


150. 


.136307 


22 


.7098103 


98.72154312222851 


64 


.98692367067588 


1 


.60 


1 


.53 


170. 


.05 


170. 


197723 


16 


.5303841 


99.33109348057584 


45 


13093652864079 


1 


.11 


1 


.16 


170. 


-06 


170. 


.188751 


17 


.0044270 


99.95114873787325 


45 


.31547803482419 


1 


.11 


1 


.16 


170. 


.07 


169. 


.485580 


17 


.7562637 


100.6056548374478 


46 


.21387973809730 


1 


.12 


1 


.17 


170. 


.08 


169. 


.477707 


18 


.0655041 


101.0199755684189 


46 


.33078105947673 


1 


.12 


1 


.17 


170. 


09 


169. 


.961853 


18 


.3975201 


101.7015952165405 


46 


.01629271848925 


1 


.12 


1 


.17 


170. 


.10 


170. 


181961 


18 


.6210308 


102.1070410738897 


45 


.89765781774504 


1 


.11 


1 


.17 


170. 


.11 


169. 


297699 


19 


.4511223 


102.8321404576463 


46 


.96444993027531 


1 


.12 


1 


.18 


170. 


12 


169. 


.478622 


19 


.8252201 


103.4316191531833 


46 


.92585236294776 


1 


.12 


1 


.18 


170. 


.13 


169 


431854 


20 


.2172241 


103.9554434974554 


47 


.09324064970371 


1 


.13 


1 


.18 


170. 


.14 


169. 


.017380 


20 


.3716507 


103.9902628410851 


47 


.51086788331664 


1 


.13 


1 


.18 


170. 


.15 


169. 


785507 


20 


.5863171 


104.6219604394106 


46 


.89467570845383 


1 


.12 


1 


.18 


170. 


-16 


170. 


062866 


21 


.1207714 


105.4821103329919 


46 


.81214457363338 


1 


.12 


1 


.18 


170. 


.17 


169. 


436630 


21 


.6631756 


105.9884280355083 


47 


.53028736542868 


1 


.13 


1 


.18 


170. 


.18 


170. 


.377396 


21 


.3885803 


105.9831545077192 


46 


.61307929154813 


1 


.12 


1 


.17 


170- 


19 


170. 


731049 


22 


.1850262 


107.2307431939022 


46 


.53162776193600 


1 


.12 


1 


.17 


170. 


.20 


169. 


.043121 


22 


.7703648 


107.4179425403126 


48 


.19269113131208 


1 


.14 


1 


.19 


170. 


.21 


169. 


.394073 


22 


.6403027 


107.3635388673966 


47 


.84569174336261 


1 


.14 


1 


.18 


170. 


.22 


170. 


449234 


22 


.9586735 


108.2095068314083 


46 


.99719888558490 


1 


.13 


1 


.17 


170. 


23 


170. 


.113739 


24 


.0381660 


109.6265273415698 


47 


.57695394709530 


1 


.13 


1 


.18 


170. 


.24 


170. 


.010727 


24 


.2990913 


109.9668119381534 


47 


.73308269329995 


1 


.14 


1 


.18 


170. 


.25 


169. 


622864 


24 


.6615162 


110.3634621605747 


48. 


.16363440273511 


1 


.14 


1 


.19 


170. 


26 


169. 


.925446 


25 


.1457653 


111.1647343607107 


48 


.01089193383915 


1 


.14 


1 


.19 


170. 


.27 


170. 


.744507 


25 


.3222599 


111.6741012150978 


47 


.33086267226461 


1 


.13 


1 


.18 


170. 


.28 


170. 


499222 


25 


.7041969 


112.1500983351095 


47 


.63318062885149 
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